ASMINT.048AUS PATENT 
HYBRID CASCADE MODEL-BASED PREDICTIVE CONTROL SYSTEM 



Reference to Related Applications 

[0001] The present application claims priority benefit of U.S. Provisional 
Application No. 60/427,143, filed November 14, 2002, titled "HYBRID CASCADE 
MODEL-BASE PREDICTIVE AND PROPORTIONAL-INTEGRAL-DERIVATIVE 
TEMPERATURE CONTROL SYSTEM FOR VERTICAL THERMAL REACTORS," the 
entire contents of which is hereby incorporated by reference. 
Field of the Invention 

[0002] This invention relates to cascade control systems having nested control 
loops for controlling a plant, such as, for example, a thermal reactor for processing of 
semiconductor substrates. 
Background of the Invention 

[0003] In semiconductor processing, the wafer size continues to increase and the 
integrated circuit feature size continues to decrease. Further increases wafer size and further 
decreases in feature size require improvements in thermal process control. The temperature at 
which wafers are processed has a first-order influence on the diffusion, deposition, and other 
thermal processes. Batch furnaces continue to play a significant role for thermal processing 
because of their large batch size and correspondingly low cost per processed wafer. A target 
in batch thermal processing is to achieve improved temperature control while maintaining 
high equipment utilization and large wafer batch sizes. The requirements of high-quality 
temperature control include a high ramp rate with good temperature uniformity during the 
ramp, fast temperature stabilization with little or no temperature overshoot, smaller steady- 
state temperature error band, shorter downtime for controller parameter tuning, etc. The 
traditional single-loop Proportional-Integral-Derivative (PDD) controllers cannot achieve the 
required temperature control performance. PID controllers with cascade or nested control 
loops have been used in attempts to provide improved temperature control. However, these 
and other previous approaches have practical drawbacks related to complexity and 
computational requirements. 
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Summary of the Invention 

[0004] The present invention solves these and other problems by providing a 
computationally-efficient hybrid cascade MBPC control system that can be implemented on a 
typical control microprocessor. In one embodiment, the hybrid cascade MBPC control system 
is a cascade-type system with nested control loops having an MBPC controller in an outer 
control loop and a conventional controller in an inner control loop. In one embodiment, the 
conventional controller is a PID controller. In one embodiment, the conventional controller is 
an H oo controller. In one embodiment, the hybrid control system uses a simplified MBPC 
control loop and a modified PID loop with a robust auto-tuner. The MBPC controller acts as 
a main or outer control loop, and the PID loop is used as a slave or inner control loop. 

[0005] In one embodiment, the hybrid cascade MBPC can be used to control a 
thermal process reactor where the MBPC controller generates the desired spike TC control 
setpoint according to both planned paddle control setpoint trajectory and the dynamic model 
related to paddle and spike TC. In the thermal process reactor, the PED loop is used to control 
the power actuator of the heater to reach the required spike control setpoint by acting as a 
local system to quickly follow changes in the spike control setpoints. 

[0006] The PID tuning parameters are relatively weakly coupled with the MBPC 
loop. The sampling time tsl in the PID control loop is preferably shorter than the sampling 
time ts2 in the MBPC control loop. In one embodiment, tsl is on the order of approximately 
1 second and tsl is on the order of approximately 4 seconds. In one embodiment, based on 
both dynamic and static models, tuning of the PID parameters is realized automatically. 
Compared with a single loop MBPC, the model order and the predictive time horizon in this 
control scheme can dramatically be reduced while the model still adequately describes and 
predicts the behavior of the actual system. In one embodiment, model derivations are done 
prior to wafer processing. 

[0007] In one embodiment, the expected temperature control range [r m j n , r ma x] is 
divided into R temperature sub-ranges flT min i, r max i], [7Wiin2> Tn^a]-- ITmin r, ^max r], where 
Tmm = Tmini , T m2LX = r max r and T m ^ x r -i = 7 m j n r ). This allows the use of linear dynamic models 
for an adequate description of the dynamic behavior within each temperature zone. 
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Corresponding to each sub-temperature range and heating zone, two kinds of dynamic linear 
models are built as: 



Pd m (t) =/nr(5/7nr(0) 
Spnr(t)=g n r(PWnr(t) 9 Pd nT (t)) 



(1) 

(2) 



where n is the heating zone number, r is the temperature sub-range, Pd is the temperature 
measured by the paddle TC, Sp is the temperature measured by the spike TC, Pw is the 
system power output, and^, r and g n ,r are linear functions. 

[0008] In one embodiment, the dynamic model of Equation (1) is used for both 
MBPC control and sofl-sensor computing, and the dynamic model of Equation (2) is used for 
soft-sensor computing and PID parameter auto-tuning of the inner control loop. 

[0009] Corresponding to the expected temperature control range, the static 
polynomial models are built as: 



where, h n are the static polynomial models. In one embodiment, the static models are used in 
limiters of the MBPC controller and for inner PID parameter auto-tuning. 

[0010] In one embodiment, the MBPC control algorithm embeds intuitive tuning 
parameters (e.g., k u , k s ) into the control law, the trajectory planner and the limiters. The 
intuitive tuning parameters can be used to improve both the dynamic control performance and 
the static control performances. The simplified MBPC control structure and fixed-time 
predictive horizon avoids the need of online matrix inversion during wafer processing. As a 
consequence, the online computing overhead is greatly reduced. In this way, the hybrid 
cascade MBPC control system algorithm can be implemented on microprocessors typically 
used in practice in the semiconductor processing industry. 

[0011] In one embodiment, a generic trajectory planer is added to the MBPC 
control loop to generate the temperature control setpoint reference trajectory. Based on the 
desired ramp rate and temperature range, the trajectory planner divides the temperature range 



Sp n =h n (Pd n ) 



(3) 
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into two sub-ranges: fast ramp; and reduced ramp. In the fast ramp sub-range, the planner 
generates the temperature control setpoints reference trajectories to enable the MBPC to 
achieve the desired ramp rate. In the reduced ramp sub-range, the planner provides at least 
one intuitive tuning parameter to control the temperature ramp speed to reach the desired 
control setpoint. Temperature stabilization time and overshoot are also controlled. This 
provides a flexible way to meet the varying temperature control requirements from the 
different processes. 

[0012] When a temperature ramp range covers more than one temperature sub- 
range, the MBPC switches its internal dynamic models so that the dynamic model, operative 
at a certain moment, corresponds with the actual temperature sub-range at that moment. In 
one embodiment, fuzzy logic switches and inference are added to the MBPC control loop to 
realize a smooth transition of the dynamic models when going from one temperature sub- 
range to an adjacent temperature range. In this case, fuzzy inference is effective to bring 
about a gradual change of one dynamic model to the other without inducing extra disturbance 
into the control system. 

[0013] In one embodiment, static limiters based on static models are embedded in 
the MBPC loop. The limiters help the MBPC to generate the correct control setpoint for the 
inner-control loop under various control cases (normal, faster/slower ramp, boat in/out, 
different load or gas flow and so on). 

[0014] In one embodiment, the PID controller in the inner control loop is 
provided with a parameter auto-tuner and/or anti-wind-up capability to enhance the 
robustness of the PID controller, and simplifies its usage. The inner PID follows the spike 
control setpoint changes that are generated by the MBPC control loop. 

[0015] In one embodiment, a software detector and control logic are included to 
detect TC measurement hardware failure. When a TC sampling failure appears, the detector 
and control logic switch on a related soft temperature sensor that is based on dynamic model 
computing. Then the soft-sensor is used to replace the real TC in its place as a control system 
input. This prevents the reactor operation from shutting down, and reduces the loss of the 
whole batch process due to the detection of one or more temperature measurement hardware 
failures. 
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Brief Description of the Drawings 

[0016] FIG. 1 shows a vertical thermal reactor with spike Thermocouples and 
paddle Thermocouples. 

[0017] FIG. 2 shows a prior art single loop paddle PID control configuration. 

[0018] FIG. 3 shows a prior art cascade PLD control configuration using nested 
PID control loops. 

[0019] FIG. 4 shows a prior art single loop MBPC control configuration. 

[0020] FIG. 5A shows a hybrid cascade control system having nested control 
loops with an outer loop and an inner loop. 

[0021] FIG. 5B shows a hybrid cascade MBPC control system having nested 
control loops with a MBPC outer loop and a conventional controller inner loop for 
controlling a thermal process plant. 

[0022] FIG. 6 shows a Hybrid MBPC control system based on a MBPC and a 
PID controller. 

[0023] FIG. 7 shows a vertical thermal reactor with the hybrid cascade 
MBPC/PID control configuration. 

[0024] FIG. 8 shows a MBPC configuration used for model identification. 

[0025] FIG. 9 shows a modeling sample data set 1 in which the spike TC 
temperatures are controlled in steps at 5 different temperature levels according to the MBPC 
configuration for model identification of FIG. 8. 

[0026] FIG. 10 shows a modeling sample data set 2, corresponding to sample data 
set 1 of Figure 9, showing the free response results of paddle TC temperatures. 

[0027] FIG. 1 1 shows a modeling sample data set 3, corresponding to sample data 
set 1 of figure 9, showing the recorded power outputs. 

[0028] FIG. 12 shows the MBPC controller structure with a trajectory planner, a 
fuzzy inference engine, and a static model limiter. 

[0029] FIG. 13 shows the PID controller structure with an auto-tuner and an anti- 
wind-up module. 
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[0030] FIG. 14 is a graph showing the temperature control setpoint generating 
principle used by the paddle TC control setpoint trajectory planner. 

[0031] FIG. 15 shows the definitions of the fuzzy sets that are used for the fuzzy 
inference engine for smoothing transition between multiple dynamic model outputs. 

[0032] FIG. 16 shows the dynamic and static relationships between the paddle TC 
and the spike TC, which illustrates the non-linearity of the vertical thermal reactor. 

[0033] FIG. 17 shows typical control results of a hybrid cascade MBPC and PID 
controller, demonstrating that in various temperature controlled ranges both ramp up 
uniformity and steady-state temperature variations are relatively small. 

[0034] FIG. 18 shows a safety control configuration 1, wherein a software signal 
based on spike TC input replaces the signal from a failing paddle TC. 

[0035] FIG. 19 is a safety control configuration 2, wherein a software signal based 
on a paddle TC input signal and power control signal replaces the signal from a failing spike 
TC. 

[0036] FIG. 20 is a graph showing the temperature uniformity during ramp-up. 
[0037] FIG. 21 is a graph showing the influence of the MBPC controller 
parameter K s for achieving different system stabilization times. 

Detailed Description 

[0038] A typical vertical thermal reactor 100 is shown in FIG 1. The vertical 
thermal reactor 100 includes a long quartz or silicon carbide process tube 110 delimiting a 
process region. A batch of wafers 152, accommodated in a wafer boat 150, placed on a 
pedestal 151 for support and thermal isolation, are inserted into the process tube 110. The 
process tube 1 10 includes an inlet 1 1 1 and an outlet 1 12 for process gas. The process tube is 
surrounded by a heating element 120 having multiple zone electric heating coils 121 to 125. 
Each zone has one or more temperature sensors. In Figure 1, each zone has a spike 
Thermocouple (TC) 130 and a "profile" or paddle Thermocouple (TC) 140. The spike TC 
produces a spike TC signal corresponding to a spike temperature. The paddle TC produces a 
paddle TC signal corresponding to a paddle temperature. The spike TCs 130 are located 
outside the process tube 110 relatively near the heating element and the paddle TCs are 
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located inside the tube 110 relatively near the wafers. The vertical reactor system 100, using 
the resistive heating element 120 to control temperature, is an inherently non-linear system 
because an electric heating element can only generate, not absorb, heat. Further, due to the 
large physical mass of the heating element 120, process tube 1 10, and wafer batch 152, and a 
correspondingly high thermal mass or heat capacity, the vertical thermal reactor 100 exhibits 
long time constants or delay times. This means that after increasing the power input of one or 
more of the heating coils 121-125, it takes a relatively long time before a new steady-state at 
a higher temperature is achieved. When the reactor temperature is higher than the control 
setpoint, cooling occurs at a rate that depends on the reactor design and its temperature, not 
on the controller. In the case of overshoot, it takes a relatively long time before the system is 
cooled down to the control setpoint again, especially at lower temperatures where the natural 
cooling rate is lower. 

[0039] Historically, the temperature in thermal reactors was controlled by a spike 
TC control loop, using a PID control algorithm. By profiling the furnace in a static mode, 
using a paddle TC, the relation between paddle TC and spike TC under static conditions was 
established and stored in a profile table. Such a profiling procedure was performed at regular 
intervals or after maintenance. Because the paddle TC gives a more relevant reading for the 
actual wafer temperature, there has been a desire to use a paddle TC control loop that would 
make the time-consuming profiling procedure unnecessary. A control configuration 200 
employing such a paddle TC control loop is shown in FIG. 2. In the configuration 200, an 
adder 210 computes an error signal Es Pd from the paddle control setpoint Pdset and the actual 
paddle temperature Pd. Based on the error signal Espa, the PID controller 220 generates a 
power output signal Pw that is provided, via a thyristor unit, not shown, to the heating 
elements 230. The tube and wafers are indicated by 240. A feedback loop 250 that includes a 
digital filter 251, provides the actual paddle TC signal to the adder 210. However, such a 
paddle TC control loop, as shown in FIG. 2, has such a strong non-linear behavior and long 
time constants that it is difficult or impossible to achieve a stable control loop with 
acceptable performance under dynamic conditions. 

[0040] A so-called cascade controller 300, as shown in FIG. 3, uses an inner spike 
temperature control loop PED controller, with the spike temperatures as inputs, to control the 
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heating elements. An additional outer loop, with the paddle TC temperatures as inputs, is 
used to generate time-dependent setpoints given to the spike temperature control loop. In the 
cascade controller 300, the spike TC temperatures are controlled by a PID controller 340 in 
.an inner loop 370. An adder 330 provides an error signal Ess p to the PED controller 340. The 
error signal Ess p is based on a spike TC control setpoint Sp se t and the actual spike 
temperature Sp. The PID controller 340 provides a power output signal Pw that is provided 
via a thyristor unit, not shown, to the heating elements 350 to heat a tube and wafers 360. A 
second PED controller 320 in an additional, outer, control loop 380 generates a spike control 
setpoint Sp set . The second controller 320 receives a paddle error signal Es Pd from the adder 
310, calculated from the paddle control setpoint Pdset and the actual paddle temperature Pd. 
Feedback loops 370 and 380 include digital filters 371 and 381 respectively to remove 
spurious data. 

[0041] A typical thermal process starts at a standby temperature at which the 
wafers are loaded into the thermal reactor. After loading, the thermal reactor heats up to the 
desired process temperature for oxidation, annealing, drive, or CVD. After performing the 
process, the thermal reactor cools to the stand-by temperature again and unloads the wafers. 
If the standby temperature, ramp up/down rate, and process temperatures are set in reasonable 
ranges, acceptable temperature control performance can be achieved during the process by 
using the cascade PID controller 300. However, optimizing the performance of the cascade 
PED controller 300 often requires significant off-line time for tuning of the controller 
parameters such as the PED parameters. Tuning of a cascade controller is often more of an art 
than a science and usually very time consuming. The best choice of tuning parameters 
depends on a variety of factors including the dynamic behavior of the controlled process, the 
controller's objectives, and the operator's understanding of the tuning procedures. For a 
cascade PID controller, the inner and outer loop tuning is strongly coupled, which adds to the 
tuning complexity. Besides the tuning of PED parameters, for dealing with long time delay, 
the outer or "profile" PID control loop of a cascade PED controller still needs a "profiling 
table" to provide constraints. Generating the profiling table involves a procedure that requires 
many hours of off-line equipment time. Off-line time cannot be used for useful wafer 
processing and is thus very expensive. 
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[0042] With the advances in modern control technology and system identification, 
more advanced control systems, such as, for example Model-Based Predictive Controllers 
(MBPC), have been developed, but these more advanced control methods are often 
computationally complex, typically requiring matrix inversion during online processing. FIG. 
4 shows a single loop control system 400 that includes a MBPC 410. The MBPC 410 
receives as input a paddle control setpoint Pdset. The controller 410 also receives as input the 
actual spike temperatures Sp and actual paddle temperatures Pd via feedback loops 440 and 
450, respectively. The controller 410 generates a power output signal Pw that is provided via 
a thyristor unit, not shown, to heating elements 420 that heat a tube and wafers 430. Feedback 
loops 440 and 450 include digital filters 441 and 451, respectively. The use of spike 
temperatures is optional. The MBPC 410 uses a complex dynamic model of the controlled 
process to compute the predictive control signals by minimizing an objective function to 
provide on-line optimization control. In the MBPC controller 410, tuning is relatively easy as 
compared to a cascade PID controller, and the MBPC 410 can compensate for delay times. In 
the MBPC 410, treatment of constraints in the system to be controlled is conceptually simple 
and multi variable control is conceptually straightforward, though computationally complex. 
The MBPC 410 typically relies on dynamic models, which are not required for the traditional 
PE) controllers. The MBPC 410 typically requires a relatively large amount of computational 
resources, particularly when constraints are considered. This can cause problems in practice 
when the control processor cannot provide the necessary computational resources. 

[0043] FIG. 5A shows a hybrid cascade controller having an outer loop controller 
501 in an outer control loop and an inner loop controller 502 in an inner control loop. The 
outer loop controller 501 provides, to a non-inverting input of an adder 530, a control 
setpoint for the inner loop controller 502. An output of the adder 530 is an error signal that is 
provided to a control input of the inner loop controller 502. A control output of the inner loop 
controller 502 is provided to a plant 509. A sensor output 505 from one or more inner loop 
sensors of the plant 509 is provided to an input of a filter 507. An output from the filter 507 is 
provided to an inverting input of the adder 530. A sensor output 506 from one or more outer 
loop sensors of the plant 509 is provided to an input of a second filter 508. An output of the 
second filter 508 is provided to an input of the outer loop controller 501. The filters 507 and 
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508 can be omitted. In one embodiment, the inner loop sensors (corresponding to the sensor 
output 505) tend to measure one or more operating parameters of a first portion 503 of the 
plant 509. In one embodiment, outer loop sensors (corresponding to the sensor output 506) 
tend to measure one or more operating parameters of a second portion 504 of the plant 509. 
In one embodiment, the first portion 503 responds to changes in the control output of the 
inner loop controller relatively more rapidly than the second portion 504. 

[0044] In one embodiment, the sensor output 505 corresponds to sensors that tend 
to respond relatively more quickly but relatively less accurately to certain desired parameters 
than the sensors corresponding to the sensor output 506. In one embodiment, the sensor 
output 506 corresponds to sensors that tend to respond relatively less quickly but relatively 
more accurately to certain desired parameters than the sensors corresponding to the sensor 
output 505. Thus, in one embodiment, the inner loop controller 502, using the sensor output 
505 is able to respond relatively quickly to certain changes in the plant 509 but relatively less 
accurately. The outer loop controller 501, using the sensor output 506 is able to respond 
relatively less quickly to certain changes in the plant 509 but relatively more accurately to 
certain desired parameters. In one embodiment, the outer loop controller 501 is configured to 
produce a setpoint for the inner loop controller 502 to improve the controlled properties of 
the plant. In one embodiment the inner loop controller 502 includes a conventional controller, 
such as a PID controller. In one embodiment, the outer loop controller includes a predictive- 
type controller. In one embodiment, the outer loop controller includes a MBPC controller. 

[0045] FIG. 5B shows a hybrid cascade thermal process plant controller 500 that 
uses a MBPC 520 in an outer control loop and a conventional controller 540 in an inner 
control loop. The cascade controller 500 is one embodiment of the hybrid cascade control 
system shown in FIG. 5A. The MBPC 520 receives as input a paddle control setpoint 
temperature Pdset and the actual paddle temperature Pd received via feedback loop 580. In 
one embodiment, the feedback loop 580 includes a digital filter 581. The MBPC 520 receives 
as input for the model calculations the actual spike temperatures via an input 572. The MBPC 
520 calculates as output a spike TC control setpoint Spset. An adder 530 calculates a spike 
error signal Es$ p by subtracting the actual spike TC temperature Sp, received via feed-back 
loop 570, provided by a digital filter 571, from the spike TC control setpoint. The spike error 
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signal Ess ? is provided to the conventional controller 540 which generates a power output 
signal Pw that is provided via a thyristor unit, not shown, to heating elements 550 to heat a 
tube and wafers 560. The sampling time tsl in the inner control loop is preferably shorter 
than the sampling time ts2 in the outer MBPC control loop. In one embodiment, tsl is in the 
range of 1 second and ts2 is in the range of 4 seconds. The conventional controller 540 can be 
implemented using a PID controller, an H controller, etc. 

[0046] As the name implies, the MBPC 520 is based on a model of the controlled 
plant. In a thermal process reactor, the MBPC 520 typically relies on several models 
corresponding to the different thermal zones of the vertical thermal reactor. The simplest 
model is a static model, describing the relation between spike TC temperature and paddle TC 
temperatures under steady-state conditions according to Equation (3). In one embodiment, the 
static model is based on 4 th order polynomial models representing the relation between spike 
TC and paddle TC temperatures over a specified temperature range. Dynamic models 
describe the dynamic response of the system. A dynamic paddle model gives the paddle TC 
temperature as a function of spike TC temperature according to Equation (1) and a dynamic 
spike model gives the spike TC temperature as a function of power output and paddle TC 
temperature, according to Equation (2). By dividing the temperature range in a plurality of 
temperature sub-ranges, a set of linear dynamic models can be obtained. This simplifies the 
required calculations. The various models are acquired experimentally from a measurement 
procedure as described below. 

[0047] FIG. 6 shows a hybrid control system 501. The control system 501 is an 
embodiment of the control system 500. In the control system 501, the conventional controller 
540 is based on a PID controller 542. 

[0048] Figure 7 shows a vertical thermal reactor system 700 where the vertical 
thermal reactor 100 is controlled by the hybrid cascade control system 501. In the system 700, 
the process tube 110 is surrounded by the heating element 120, comprising multiple zone 
electric heating coils. Each zone has the spike TC 130 and the "profile" or paddle TC 140. 
The spike TC is located outside the process tube relatively near the heating element and the 
paddle TC is located inside the tube relatively near the wafers. A paddle control setpoint Pd se{ 
and the actual paddle temperatures Pd are provided to an MBPC controller 720 
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(corresponding to an embodiment of the MBPC controller 520), which generates a spike 
control setpoint Sp set . An adder 730 computes a spike TC error signal using the spike control 
setpoint Sp se t and the actual spike temperatures, provided to the adder 730 via an inverter 
732. A PID controller 740 uses the spike error signal to generate a power output signal that is 
provided to a power actuator 750 to provide power to control the heating element 120. 
1 . Experiment Design and Data Acquisition for Model Identification 

[0049] The identification test design plays an important role in a successful model 
identification and MBPC design. Current practices of identification methods are to use single 
variable step or finite impulse tests for MBPC model identification. The tests are carried out 
manually. The advantages of these methods are that the system dynamic responses are 
described in an intuitive manner. One drawback with the step or finite impulse tests is that 
the data from these tests may not contain enough information about the multivariate 
characteristics of the process because the step or pulse signals may not induce enough 
dynamic behavior of the process. A second drawback of these step or impulse tests is that 
they can be very time consuming. To avoid the above-mentioned drawbacks of common 
identification methods, during the identification procedure for the MBPC 720, the inner-loop 
PID control loop 740 is actively used. The PED constants used during the model identification 
are based on previous control experiences of the vertical thermal reactor. Then, using this 
inner-loop PED controller 740, the system identification can be carried out automatically by 
using a model identification and data acquisition recipe. The signals inducing dynamic 
behavior of the system are real control signals, and the conditions are similar to real process 
conditions. In this case, the PED controller 740 helps to keep the spike TC temperatures 
within their limits. The models based on these inner closed-loop data enhance the 
performance and stability margins of the system 700. 

[0050] In one embodiment, the modeling data acquisition is achieved by using the 
control scheme shown in FIG. 8. In FIG. 8, a temperature control setpoint for model 
identification is provided to a trajectory planning unit 800, which creates a time-dependent 
spike TC control setpoint temperature Sp se t such as, for example, a controlled ramp-up rate. 
An adder 810 calculates a spike error signal Es, using the spike control setpoint Sp set and the 
actual spike temperatures, received via a feedback loop 850, provided with a digital filter 
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851. A PID controller 820 uses the spike error signal Es to calculate a power output P\v 
signal that is provided via a thyristor unit, not shown, into the heating elements 830 to heat a 
tube and wafers 840. 

[0051] An example identification control process sequence (i.e., a "recipe") for a 
vertical thermal reactor starts at room temperature and ramps up at a ramp rate (that is, a time 
rate of change) of 10 °C/min, stabilizing the spike temperature for approximately 45 minutes 
at 200, 400, 600, 800 and finally 1000 °C, using a PID controller in the configuration of 
FIG.8. In order to prevent slip, the ramp rate during the last step is 5°C/min. During 
execution of the model identification and data acquisition recipe, the resulting actual spike 
and paddle temperature signals (shown in FIGS. 9-11, respectively) and the power output 
signals are recorded and stored for modeling. 

[0052] The data collected from the recipe is divided into five data subsets, 
corresponding to the five stabilization temperatures of the model identification and data 
acquisition recipe. Each subset starts at the beginning of a ramp up period and ends just 
before the beginning of the next ramp up period. 

[0053] In contrast to prior art model identification methods where data is acquired 
in open loop control, in the present invention the data is acquired in closed loop control as 
shown in Figs. 7 and 8. Therefore, the collected data is typically free of spurious data, so little 
or no data pre-processing is required to remove spurious data points. 

2. Model Structure 

2.1 Static Models 

[0054] For each zone, a static model is derived. For each stabilization temperature 
in the identification recipe, at least one value is extracted both for the spike TC temperature 
and the paddle TC temperature. These pairs of values, 5 in the example shown in FIGS. 9-11, 
are used to estimate the parameters for the static model using a parameter estimation 
technique, such as, for example, a polynomial fit, least-squares fit, etc. In one embodiment, 
the model equation for the static model is: 



-13- 



= *»o +s mi Pd H +s n2 Pd 2 n +... + s n{q . l) Pdr i +s nv Pd< n 



(4) 



where n is the zone number, s nq are the static models parameters to be determined, and q is 
the order of the static model. The model according to equation (4) gives an adequate 
description of the relationship between the spike temperature Sp and the paddle temperature 
Pd over the desired temperature range. 
2.2 Dynamic Models 

[0055] For each zone, a dynamic model is derived from each data subset for the 
spike TC and the paddle TC. This means that in the case of five zones and five temperature 
sub-ranges, according to the present example, 25 dynamic paddle linear models and 25 
dynamic spike linear models are used. It will be clear to one of ordinary skill in the art that 
any number of thermal zones and temperature sub-ranges can be selected, depending on the 
circumstances. In one embodiment, a linear least-squares algorithm is used. The model 
equations for the dynamic linear models are: 

Pd nr (t) = f nr (S Pnr (t)) 

= a nrl Pd nr (t-l) + a nr2 Pd nr (t-2) + ... + a nrl Pd nr (t-l) (5) 
+ Ko + bnnSPn, (t - 1 ) + b nr2 S Pnr (t - 2) . . . + b^Sp™ (t - m) + e lu _ 



= c nrX S Pnr (t - 1) + c nr2 Sp nr (/ - 2) + . . . + c m Sp„ (t - x) 

+ d nr0 + d nrX Pw nr (t - 1) + d nr2 Pw nr (/ - 2) + . . . + d nrx Pw nr (t - y) 

+ PnrxPd nr (* - 1) + P n nPd nr (t - 2) + . . . + Pm Pd„ (t - Z) + e„ rs 



where n is the zone number, r is the temperature range number, / is the discrete time index, 
a nr i and b nrm are paddle model parameters, / and m are model orders of the paddle models, 
and c nrx , e?nry and p nrz are spike model parameters, Further, x, y, and z are model orders of the 
spike models, and e nrp and e nK are model errors or disturbances. Typically, a first or second 
order approximation {m = 1 or 2) is adequate for the model of equation (5) whereas the order 
for the model of equation (6) is typically 20 or more (e.g., / = 28) for an adequate description. 
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[0056] Model validation can be provided by visual comparison of the measured 
and calculated model output, simulation, residual analysis, cross-correlation analysis, etc. 
3, Mathematics of Model Extraction 

[0057] The methods used for parameter estimation for both the static models and 
the dynamic models involve solving a linear least-squares problem (LLS). This can be done 
via the so-called normal equations, or via a QR-decomposition. The method via QR- 
decomposition typically requires more calculations, but tend to be numerically more stable. 
For most cases, solving the normal equations gives good results, but for higher orders it is 
safer to use the QR-decomposition. 

[0058] The method of solving the LLS problem is described below. Given a 
system of equations defined as A x = b , where A is a matrix and jc and b are vectors, the 
linear least-squares problem is to find a vector x that minimizes 

p(x)=\\A-x-b\\ (7) 

If the matrix A is non-singular (i.e., if the inverse of A exists), this problem has a unique 
solution. 

3.1 Solving LLS via the Normal Equations 

[0059] The direct method of solving a LLS problem is via the normal equations. 
If a function has a minimum, its derivative at that minimum is zero. Thus (p(x) has a 
minimum where d cp(x)/dx = 0 or: 

A T A-x-A T -b = 0 (8) 
where A T is the inverse matrix of ,4. This results in the so-called Gaussian normal equations: 

A T A x = A T b (9) 
From these equations, x can be calculated as 
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x = (A T AYA T b (10) 
3.2 Solving LLS via OR-decomposition 

[0060] Another method of solving a LLS problem is via QR-decomposition. With 
this algorithm the matrix A is expressed as the product of an orthogonal matrix Q and an 
upper triangular matrix R. 



A = QR 



01) 



[0061] The QR-decomposition of the matrix A can be computed by calculating 
the Householder reflection H for each column of A. The Householder reflection Hu of the k th 
column can be calculated as: 



(12) 



™k = 



a 



a 



(13) 



H k =I-2wy k 



(14) 



where a-^ is the element of matrix A at the I th row and / h column, and n is the number of 
columns and I is the unity matrix. 

[0062] The matrices Q and R can be calculated as 



Q = H,H 2 //„_, 



(15) 
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R = H n _ } H n _ 2 H t A 



(16) 



For a matrix A of m rows and n columns, where m is greater than n, the last m-n rows of the 
upper-triangular matrix R are completely zero. 

[0063] Once the matrix A is expressed as A=QR, the LLS problem can be written 

as: 



\Qy = b 

[Rx = y 



(17) 



Q is orthogonal, so 

y = Q T b 



(18) 



[0064] With y and R known, x can be calculated. Since R is upper-triangular, this 
is simply done via backward substitution. By way of example, the solution is shown for a 3 
by 3 matrix in the equation below. 





r 12 








"y." 
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x 2 




y 2 


. r 31 
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0_ 


. X 3_ 




_y 3 _ 



(19) 



This system can be solved via backward substitution starting at the bottom row. 



r 3i -*i = J>3 



r 2i +r n -x 2 = y 2 =*x 2 = yjr 22 -r 2x x, = y 2 /r 22 -r 2x ■ yjr 3l 
r„ +r 12 -* 2 +r 13 -* 3 => 

x 3 = yJ r n ~ r w '*> -r n ■x 2 =yjr n -r u -yjr^ -r n -(y 2 /r 22 -r 2x -yjr ix ) 



(20) 
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[0065] When using the QR-decomposition for solving a LLS problem, the 
orthogonal matrix Q is usually not explicitly calculated. Instead, R and y are calculated in a 
recursive algorithm, with initial conditions 



[0066] Next for each column of R, the vector w is calculated. With this vector, R 
and y are updated according to the following formulas. 



where Rk is the k column of R. This is repeated for all columns of R. The solution of the 
LLS problem can then be calculated via backward substitution as described before. 
4. Extraction of the Dynamic Models 

[0067] The dynamic models (shown in Equations (5) and (6)) for the MBPC can 
be represented by the following equations: 



(22) 



y 




nrp 



(23) 



(24) 



where Pd mr and Sp nr are model predictive outputs, and 
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1 

Sp nr (t-V 
Sp nr (t-m) 



and 



0 = 



itW 



'nrO 



*U<)= 



1 

-1) 



Pd nr (t-z) 



and 



0... = 



'nrl 



[0068] With the model structure defined above, given the model orders and a set 
of input and output data, the parameters of the model are found by minimizing a so-called 
loss function. An often-used loss function is the summed squared error: 



(25) 



where N is the number of input and output samples and e is the prediction error vector, 
defined as the difference between measured and predicted output: 



e pd (t\e) = Pd nr (0 - Pd nr (t\0) = Pd ttr (0 - f nrp {t)9 
e sp (t\e)=S Pnr (t)-Sp nr (t\0) = s Pnr (0-C W„, 



(26) 
(27) 
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For simplification, the matrix O and vector Y are used, defined as: 





'«,(!)" 














0 = 


t(2) 


or = 


C(2) 


y = 


«t(2) 


or = 


S Pnr {2) 




AW 















(28) 



The loss function now becomes: 



j N {e)=\Y-®o\ 



(29) 



where 9 is equal to ^ or # nrs . 

[0069] Minimizing this loss function can be accomplished by solving a linear 
least-squares problem. 
5. Extraction of the Static Models 

[0070] The parameters of the static models shown in Equation (4) are obtained by 
polynomial fitting using groups of input and output data. Thus, Equation (4) can be rewritten 
as: 



S Pn (*) = s n0 + s n] Pd„ (k)+ s n2 Pd n (kf + ... + s n Pd„ {k)" 



(30) 



where k is the k th value of the input and output sequence. 

[0071] The identification problem can now be formulated as follows. Given input 
and output signals Sp n =[Sp n (l), Sp n (2), ... , Sp n (k)] T , Pd n =[Pd n (l), Pd n (2), Pd„(k)) r , and 
model order q, find appropriate values for parameters S n =[s nq , s n i, s n o] T First a 
Vandermonde matrix V is constructed: 



V = 



PdM - M„(l) I 
Pd„(k) 9 - PdM 1 



(31) 
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Next, the parameters are estimated by solving the following least-squares problem. 









- Pd„(\) 1 




_SpM 




r*M ■ 







which can be written as 

Sp„=VS„ (33) 

Here the identification problem has become a linear least-squares problem, thus the 
parameters of a polynomial fit can be calculated as 

S„=(v T vYv T Sp„ (34) 

Alternatively, the QR decomposition algorithm can also be used to find the parameters for the 
polynomial fit. 

[0072] Model validation can be provided by visual comparison of the measured 
and calculated model output, simulation, residual analysis, cross-correlation analysis, etc. 
6. MBPC Structure 

[0073] The internal structure of one embodiment of an MBPC corresponding to 
the MBPC 520 from FIG. 5 is shown in more detail in FIG. 12 as an MBPC 1200. The 
MBPC 1200 includes an MBPC algorithm module 1230, a trajectory planning module 1220, 
an MBPC outputs fuzzy inference module 1240, and an MBPC static model limiter 1250. 
The MBPC algorithm module 1230 includes a modeling module 1231, which performs the 
actual modeling based on the dynamic model of Equation (5), and an optimizer module 1232. 

[0074] Inputs to the MBPC controller 1200 are the paddle control setpoint 
temperature Pd set and the actual paddle temperatures P d . The paddle control setpoint 
temperature is provided to the trajectory planning module 1220 and the actual paddle 
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temperatures P d are provided to a memory 1210 for storing past inputs and outputs. The 
memory 1210 provides input to the MBPC algorithm module 1230. Additional input for the 
models include actual spike temperatures. The Trajectory planning module 1220 generates N 
paddle control setpoints Pd set (l...A0 distributed over a predictive horizon, where Pd set (l) is 
the control setpoint for the present moment and Pd set (A0 is the most future predicted control 
setpoint. These control setpoints Pd se t(l...A0 are provided to a first input of an adder 1220. 
Further, the modeled paddle values Pd fr (\...N) , which are provided as output by the MBPC 

control algorithm module 1230, are provided to a second input of the adder 1222 via a line 
1233. The adder 1220 calculates error signals Es(l...A0 which are provided to the Optimizer 
module 1232 of the MBPC algorithm module 1230. The optimizer module 1232 optimizes 
the model output by minimizing a cost function 1235 as represented by equation (35), using 
constraints 1236. The least-squares error between the modeled predicted paddle control 

setpoint temperatures Pd /r (l...N) and the actual paddle control setpoint temperatures 

Pd set (l...N) from the trajectory planner 1220 is minimized over the predictive horizon. The 
predicted paddle control setpoint temperatures are optimized by using the disturbance model 
(the last term in equation (35)) so that the predictive values approach the actual values. 

[0075] The spike correction value ASp is calculated, according to equation (45). 

The modeled values Pcl fr (l...N) are provided to the memory 1210. The spike correction 

value ASp is provided from the MBPC algorithm into a spike output calculation module 1212 
to calculate the modeled spike control setpoint Sp sei (l) according to equation (46). The 
modeled spike control setpoint Sp se t(l) is provided to the MBPC outputs fuzzy inference 
module 1240. The modeled spike control setpoint value is provided to the Output limiter 
1250 where the output is limited according to equation (54). The algorithms will be discussed 
in further detail below. 
6.1 MBPC Algorithm 

[0076] Based on the dynamic linear models described in Equation (5), the 
predictive control algorithm calculates the control strategy ASp s a(0 by minimizing the cost 
function J, defined as: 
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/ = (t + k)- Pd(t + k | t)f + Y,k s [*Sp{t + k- 1)] 2 (35) 

k = ] k = \ 



where N and N u are the prediction horizon, k u and k s are weight parameters, and Pd SQi (t + k) is 
the k ih paddle control setpoint generated by the trajectory planner. Further, Pd(t + k\t) is the 
k {h model predictive output at time t, which can be considered as the combination result of 
two separate contributions: 

Pd(t + k\t) = Pd f ,(t + k\t) + Pd fo (t + k\t) (36) 

where Pd fr (f + k \ t) is the free response, and Pd fo (t + k\ t) is the forced response. Among 
them, Pd fr (t + k \ t) can be computed directly from Equation (6) as: 

/ m n 

Pd fr (t + k |0 = £ ~i* + >t|0 +2 - y + * 1 0 + Yj d r e( J -r + k\t) (37) 

1=1 7=0 r=l 

n 

where ^</ r e(* - r + A: | /) is the disturbance model output, and 

e(0 = [«(#) - Pd fr (/)] - « - 1) - $/>(<)] • 
Then, P</ /<( (f + Ar 1 1) can be calculated as: 

Pd f0 (t + k\t)= Xg^ASpit + i\t) (38) 

1=0 

where g\ is the coefficient of the module step response of the model from Equation (5), 
which can be obtained as: 
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1-1 



S.=Z fl y*.-; + I>y =0 when < 39 > 



where Oj and b } are the coefficients from Equation (5). 

[0077] By using matrix notation, Equations (38), (36) and (35) can be rewritten 

as: 

Pd f0 =GASp (40) 

Pd = Pd fr+ Pd f0 
= Pd fr +GASp 

J = k u (Pd SH -Pd) T (Pd set -Pd) + k s ASp T ASp 

= k u [(Pd sel -Pd fr )-GASp] T [(Pd sel -Pd fr )-GASp] + k s ASp T ASp 

dJ 

[0078] Minimizing J with respect to ASp, i.e., = 0, then ASp gives the 

d(ASp) 

optimal solution: 

ASp = (k u G T G + k s iy i G T k u (Pd se ,-Pd fi ) (43) 

where I is a identity matrix, and the matrix (k„G T G + k s I) to be inverted has N u x N u 
dimension. The value of N u need not be 1, but it is instructive to note that when N u = 1, then 
Eq. (43) is simplified as the scalar control law: 



k„ Z g k [Pd se , (t + k\t)-Pd, r (t + k\ 0] 
ASp(t\t) = —^ (44) 
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Substitute Eq. (39) into (44), then: 



*.Z 



* ( 



I OjSi-1 + I f Pd s« (t + k\t)-Pd„(t + k\t)] 



k k-\ \ 



ASp(t | 0 = 



(45) 



N ( k k-\ V 



[0079] In one embodiment, the MBPC employs a "receding horizon" control 
principle, where only the first element ASp(/|/) is required to compute the MBPC output: 



At the next sampling instant the whole procedure is repeated. 

6.2 Trajectory Planning 

[0080] To control the desired speed to approach the paddle TC control setpoint, a 
trajectory planer is used. In one embodiment, the trajectory planner reduces the ramp rate 
when the actual paddle temperature Pd approaches the desired paddle control setpoint Pd w to 
within a range d: 



where r p is the desired ramp rate, £is the stabilization range and a is the trajectory reference 
time constant. An example trajectory reference is shown in Figure 14. Both 8 and a are 
related to, and depend on, the ramp rate r p , and are separately defined as: 



Sp set (') = Sp set (t-l) + ASp(t\t) 



(46) 





(47) 




^[S min ,d max ] and k r <1 



(48) 
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where k r is a proportional constant, and k s is a constant for the stabilization time control that 
is also used in the control law in Equation (44). 

[0081] The parameters 3 and a can be selected and adjusted to provide optimum 
control under a variety of circumstances by setting the values of the constants k r and A; s . It can 
be seen from FIG. 14 that by tuning k s , different trajectories are achieved, which provides a 
natural way to control the aggressiveness of the controller. Increasing k s is analogous to 
increasing a, which leads to a slower but more robust controller. Thus, k s is a more direct and 
intuitive tuning parameter for system stabilization than a. Similarly, k r is a more direct and 
intuitive tuning parameter for system stabilization than S. 

6.3 MBPC Outputs Inference by Using Fuzzy Logic 

[0082] During controlling temperature ramp up/down to a desired temperature 
level, the temperatures may cross several different temperature sub-ranges, particularly, when 
the temperature is close to the defined edges of a temperature sub-range. The MBPC needs to 
switch from the dynamic models that are valid for one temperature sub-range to the dynamic 
models that are valid for another, adjacent, temperature sub-range. In one embodiment, a 
fuzzy logic software algorithm is used to ensure that the MBPC controller has a smooth 
transition of its outputs, without introducing disturbances into the control system when it 
needs to switch between the models. 

[0083] The fuzzy sets are defined as shown in FIG. 15. According to 
Takagi/Sugeno fuzzy rule format, the fuzzy rules are written as: 

R\: IF paddle TC is in temperature sub-range S m i n Then Sp 1 set (0 = Sp r \ (t) 
ALSO 

R2: IF paddle TC is in temperature sub-range S2 Then Sp 2 set (0 = Spa (t) 
ALSO 
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ALSO 

R m : IF paddle TC is in temperature sub-range S max Then Sp m se t (0 = S/?rm (0 

where, S m i n , S2, . . S max are fuzzy sets defined in FIG. 15, Sp™ are the crisp values of MBPC 
outputs Equation (46) from the different models. 

[0084] By using fuzzy inference and a center of area (CO A) defuzzification 
operator, the outputs of the MBPC are calculated as: 



5>i 

*:w aJ h 

!>.(*.) 



(50) 



where « is the zone number, m is the fuzzy rule number, Sp" m (f ) is the spike control setpoint 

for zone n, /j\ (i=l, 2, m) are the fuzzy membership functions that are defined in FIG. 10 
and can be mathematically described as: 



ft (*«.) = 



1 



2AT 



Sp sel e(0, 7V-A7-) 
Sp^^T^-AT, T min +AT] 
Sp^eiT^+AT, +co) 



(51) 



where [Jl™, ?max] is the desired temperature control range, AT is a parameter that 
characterize the membership functions. 



M 2 =... = Hj(Sp set) ) = 



0 

S Psel -(T H -AT) 



2AT 



(Tj + AT) - Sp s 



0 



2AT 



Sp sel e(0, Tj_ x - AT) 

Sp m e[T„-*T, Tj. l+ AT) 
S Pstl e[T M +AT, Tj-AT) 

S Psel e[Tj-AT, T j+ AT] 
Sp sel e(Tj+AT, +co) 



(52) 
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where j = 2, . . m-1 are the sub-temperature control ranges , and Ti = T m i, 



M m (Sp set ) = < 



0 



Sp set e[T n „ x -AT, T max +AT) 
Sp^eiT^+AT, +co) 



(53) 



2AT 

1 



6.4 MBPC Output Limiters 

[0085] The hybrid cascade MBPC and PED control scheme has an important 
advantage over single loop systems in that it reduces the disturbance effects of inner-loop 
process variability. However, such a scheme also presents some conceptual difficulties, and 
can result in poorer dynamic performance if the inner loop PID is not fast enough to follow 
the control setpoint generated by the MBPC. Moreover, model mismatches can appear in 
cases where the reactor system is changed or disturbed (e.g., boat in/out, when the door of the 
reactor is open, different ramp rates and load sizes, etc.). In these cases, the conceptual 
difficulties arise from the fact that the MBPC needs to know its control effort limits in order 
to function properly. To deal with these difficulties, limiters based on static models are added 
to the MBPC control loop. The limiters are defined as: 



where Sp n is computed by Equation (30), AT is an adjustable temperature constant for 
compensation of the model outputs mismatches, and k u is a tuning parameter that is also used 
in the control law Equation (45). Tuning k u , can improve the temperature uniformity (shown 
in FIG 20). 

[0086] It can be seen from FIG. 16 that both dynamic and static relationships 
between the paddle TCs and the spike TCs are changed with the different temperature levels. 
The dynamic relationships are different from the static relationships, especially for lower 



SpZ" 

split) =\sp';At) 



Split) < Sp™ =-\Sp„it) + k u AT\ 

- 1 Sp n it) + k u AT | < Spl )<\Sp„(t) + k u M\ 

Split) >Sp™=\Sp a (t) + k u AT\ 



(54) 



-28- 



temperature ranges. These differences result from the temperature non-linearity of the 
thermal reactor. 

[0087] By using Equation (54), some uncertain factors caused by the system non- 
linearity and disturbances are limited in the MBPC control loop, which ensures that the 
MBPC can always generate the reasonable control setpoints for the inner PID control loop, 
enhances the stability, control margins and robustness of the hybrid cascade MBPC control 
system. 

7. PEP design 

[0088] The PID loop is an inner control loop that works in concert with the outer 
MBPC control loop. Tuning parameters (e.g., control gains) of the PID controller are based 
on the accumulated control experiences and the open-loop identification analysis of the 
vertical thermal reactor. The PID controller is used for both modeling (FIG. 8) and normal 
control (FIG. 6, 18 19). Its structure is shown in FIG. 13. The PID algorithm is 
mathematically represented by: 

p w (t) = if E s (t) + Iout(t) + Dout(t) jxG 
S 

= l f E s O)^ifE s (t) + ^-D s (0J J+ ^-D s (t)}xG 
= k p (t){ E s (t) + k, (f)£ / E s (t) + k d (t) D s (t) J j + k d (t) D s (t)}xG 
(55) 
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where : 




Fw(t) = 


- power output 


d = 


= proportional band 


E s {t) 


= op 5et (t)-bp(t) = spike error 


= 


*i < • < 
= spike setpoint 


Sp(t) -- 


= spike sampling temperature 


rrt 

T s 


= sampling period (1 second) 


T, 


= integral time constant 


T d 


= differential time constant 


DM) - 


= Sp(t -1) - Sp(t) = spike TC change rate 


G 


= General gain for sepecial power control 



The first, second and third terms on the right-hand side of Equation (55) represent the 
proportional, integrative and differentiating actions respectively. 

[0089] The structure of a PID controller 1300 is shown in Figure 13 and reflects 
the structure of Equation (55). In the PED controller 1300, an error signal Es is provided to a 
limiter 1302, mathematically represented by equation (57) below. After passing through the 
limiter 1302, the limited error signal is provided to a differentiating action module 1320, to 
an integrating action module 1340, and to a proportional action module 1360. In the 
differentiating module 1320, the rate of error signal change is determined in block 1322. 
Then the differentiating action is calculated in block 1324 using a k d value, which is variable 
and calculated in block 1328. The calculated differentiating action passes the output limiter 
1326, mathematically represented by Equation (59). 

[0090] In the integrating action module 1340, an integration constant kj is 
calculated in block 1348 and applied in block 1342. A summation is carried out in block 
1344. Then the calculated integrating action passes through the output limiter 1346, 
mathematically represented by Equation (58). In one embodiment, the output of the 
differentiating action module, Dout, is also used as an input for calculating the integrating 
action, as shown by feedback line 1331. The purpose of this feedback is to achieve improved 
control during ramp-up and ramp-down. The outputs of both the differentiating action 
module 1320 and the integrating action module 1340 are provided to the proportional action 
module 1360. In the proportional action module 1360, the proportionality action is calculated 
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in block 1362, using a proportionality constant k p calculated in block 1366. Then the output 
signal Pw passes through output limiter 1364, represented by Equation (60). 

[0091] The PDD parameters k d , k\ and k p are calculated in blocks 1328, 1348 and 
1368 according to the formulas: 



MO 
MO 



_ p 



r p (0T ds (0 



= K max"*! »x(L,-*«,(0) 



T^ZSp sel (t)ZT m 



(56) 



where: 



A: 

/>_min 



it 

/_min 



mm 



7 1 -7 



max mm 



**.--(^o+*/.x( ^ (7 '- h) -^)) 

■* mi 



mm 



*■ max * mm 



where Sp SQt (T) is calculated by using the static model according to Equation (3), and k p0 , £ p i, 
kp2, k\o, k{\ 9 and k\ 2 are predetermined constants relating to system gain and time constants. 
Tds is a delay time. T m \ n and T m2LX are the lower and upper temperature boundaries of the 
temperature control range. After dynamic and static models are identified, the k p ma x, & p _ m j n , 
£i_max, £i_min and Tds can be re-determined and modified by performing an analysis of both the 
dynamic and static models. For a given thermal reactor, these parameters can be 
predetermined in the design and development phase without additional on-line tuning for 
individual reactors after manufacturing or during use. 

[0092] For anti-wind-up and integral saturation, the following linear limiters are 
included in the PED controller 1300, and defined as: 
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EM) = 



Es_min Es(t) <Es_min 

Es(t) Esjnin < Es(t) < Es_max 

Esjmax Es(t) > Es_max 



(57) 



Iout(t) = 



Ioutmin 
Iout(t) 
lout max 



Iout(t) <Ioutjniin 

Ioutmin < Iout(t) < loutjtnax 

Iout(t) > lout max 



(58) 



Dout(t) = 



Doutjmin Dout(t) <Dout_min 

Dout(t) Doutjnin < Dout(t) < Doutjtnax 

D out max Dout(t) > Dout max 



(59) 



Pw(t) = 



Pwmin Pw(t) < pwjmin 

Pw(t) Pw min < Pw(t) < Pw max 

Pwmax Pw(t) > Pw max 



(60) 



[0093] By using Equations (57)-(60), the dynamic response of the inner-loop PED 
is stable and provides the desired speed of response. The control results can be seen from 
modeling data sets in FIG. 9 and FIG. 10. The design of the outer loop MBPC based on the 
dynamic models derived from the closed-loop data controlled by the PID controller greatly 
increases the temperature system control margins and stability. 
8. TC Measurement Failure Detection and Remedy 

[0094] It can be seen from FIG. 1 that the control system has two groups of 
temperature measurement sensors: 1) paddle TC; and 2) spike TC. In one embodiment, each 
group includes five sampling circuits, each circuit corresponding to a specific heating zone. 
During the batch process, the software detectors in the temperature control system monitor 
the TC measurements and determine whether there are failures or not. If one or more TC 
hardware failures are detected, then the control logic switches the TC input(s) from real 
sampling measurements to corresponding soft-computing values that are calculated based on 
dynamic models in Equations (5) or (6). Using the soft-sensors, the soft-sensor temperature 
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control quality is substantially equivalent to the normal temperature control, which prevents 
the reactor operation from shutdown due to TC measurement hardware failures. The process 
under the procedure can continue without the loss of the whole batch process. 

[0095] The MBPC safety control configuration for paddle TC failure is shown in 
FIG. 18 which is similar to FIG. 6, with similar parts indicated with similar reference 
numerals. Additionally, FIG. 18 shows a soft-sensor paddle 1810 that models paddle 
temperatures based on spike temperatures. Further a TC failure detection module and switch 
1 820 are provided. When a failure of a paddle TC is detected, the actual paddle temperature 
for that paddle TC is no longer provided as input to the adder 510, but the modeled paddle 
temperature, based on the spike values, is used. 

[0096] Similarly, the MBPC safety control configuration for spike TC failure is 
shown in FIG. 19. FIG 19 shows a Soft sensor spike 1910, which models a spike temperature 
based on paddle temperature Pd and power output Pw. Further, a TC failure detection and 
switch 1920 is provided. When a failure is detected for a spike TC, the actual spike 
temperature for that TC is no longer provided as input to the adder 530 and to the MBPC 
controller 520, but instead the modeled spike temperature is used for that purpose. 
9. MBPC Parameter Tuning 

[0097] The presence of the inner-loop PUD control with outer-loop MPBC control 
makes the control system robust such that in many cases parameters tuning only needs to be 
performed during the design phase. For small changes in the furnace system, additional 
parameter tuning can be omitted in many cases. Since both the dynamic and static models 
are derived from the closed-loop data controlled by the inner-loop PID controller, slowing 
down or speeding up the inner-loop does not seriously degrade the performance of the outer 
MBPC loop. In this case, the inner and outer loop tuning are not strongly coupled. In the 
design of the inner loop, parameter tuning can be realized automatically (shown in Equation 
(56)). For control flexibility, in the outer MBPC loop design, two parameters (k u and k s ) are 
provided for the adjustment of the dynamic ramp-up temperature uniformity (shown in FIG. 
20) and the stabilization time (shown in FIG. 21). Both k u and k s are embedded into the 
control law expressed by Equation (45). Further, k u is also included in the MBPC output 
limiters Equation (54) and k s in the trajectory planner Equation (49). Both k u and k s behave 
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like time constants. Thus increasing/decreasing k n or k s is analogous to changing the control 
system time constants, which causes the system response to slow down or to speed up. The 
main effect of k u is on dynamic temperature uniformity during ramp-up or ramp-down, and 
the main effect of k s is on stabilization time control. The tuning rules can be summarized 
and simplified as: 

• Increase k u , the ramp up rate slows down; 

• Decrease k Uj the ramp up rate speeds up; 

• Increase k s , the stabilization time becomes longer; 

• Decrease k s , the stabilization time becomes shorter. 

[0098] Although the present invention has been described with reference to a 
specific embodiment, other embodiments occur to those skilled in the art. For example, the 
hybrid cascade MBPC can be used to control many linear and/or non-linear plants, not just 
the vertical thermal process reactor. It is to be understood therefore, that the invention herein 
encompasses all such embodiments that do not depart from the spirit and scope of the 
invention as defined in the appended claims. 
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